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This paper presents a method of solving for Fourier coefficients where the dependent variable 
can be expressed as a piecewise continuous function, when various conditions of continuity and 
smoothing are assumed. An example is included to show the effect of smoothing in the region of a 
discontinuity for a system composed of two materials that exhibit a discontinuity at their interface 
and surrounded by a third material which does not have a discontinuity. An advantage to be gained 
from smoothing is an increase in the convergence of a finite Fourier series representation of a piece- 
wise continuous function in the region of the discontinuities. 
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1. Introduction 

Many time-independent problems of a practical nature involve the solution to the partial 
differential equation 

d 2 v d 2 v d 2 v__~ 

dx 2 dT dz 2 

in rectangular coordinates, or 

d 2 v 1 dv 1 d 2 v d~v _ , 

dr 2 r dr r 2 dcfr dz 2 

in cylindrical coordinates. From the general solution, one particular solution is in the form 

v = B + Cz + y An sin— , (2) 



where the Fourier coefficients, A n = A n {x, y) for rectangular coordinates, or A n = A n (r, (p) for cylin- 
drical coordinates, are determined from boundary conditions in the region ^ z ^ co, and fixed 
coordinate positions for x and y, or r and <£>, with the condition that other necessary solutions of 
(1) disappear on this boundary. Also, B + Cz must satisfy v at the end points, z = and z = to. 

From an experimental point of view, the dependent variable, v, can be approximated from 
measured values at several z positions on the boundary by a relationship, v = g{z). The Fourier 
coefficients are usually simply determined if the function g(z) and its derivatives are continuous 
over the region ^ z ^ co. For various applications, the variable v must be expressed as piecewise 
continuous functions in this region, with apparent discontinuities at discrete points. If it is known 
that the discontinuities do not exist, then an appropriate smoothing function can be employed. 

The purpose of this paper is to solve for Fourier coefficients, when it is assumed that the 
dependent variable v can be expressed by polynomials representing discrete portions of the region 
^ z ^ oj, and that conditions of continuity and smoothing exist. An application of functions 
developed to a hypothetical problem is also presented. 
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2. Analysis 

On a fixed surface lying in the region ^ z ^ co, the dependent variable v is defined 

v = gt(z); zi-i <z< zi, 

where i— 1, 2, 3, . . . m, Zo = 0, and z m = (o. The functions gi(z) and their derivatives may 
have discontinuities at the points Zi, i = 1, 2, . . . m — 1. Substituting in (2), multiplying through 
by sin A;77z/co, and integrating between the specified limits gives the following equations 

i 1 rC / TTT I i T1TT7 rCTT? 1 

[gi(z)—B — Cz] sin dz = \ A n sin sin dz. 

J*t-l W n J*/-i W W 

Adding the above equations 

f w 717TZ ^ 77 - z f w A:7TZ m C z i klTZ 

\ A n sin sin- dz — —\ (B + Cz) sin cfeH-V giU) sin dz, 

~ Jo w oj Jo co £lJZi-l w 

from which 

2 9 m f Z i T17T7 

A n = — {Cfti(-l)»-B[l-(-l)»]}+-y #(*) sin— c/z. 

^vr co^-jj^j to 

Repeated integration by parts yields 



. tt77"Z 

sin 



[ ( , . TITTZ tof ,, . gj (z) g?(z) I 

#(z) sin dz = ^ — \riftW~7 — \2~ + 7 — \J~ ' * -\ 

J co (]UL\ I [ZL?L| [ML] J 

\ Co/ \ ft)/ \ ft)/ 



nirz 

cos ~ r M gru) , #(*) 

l#w — 7 — T2 + 7 — \1 

M7T L (niry fniry 

co \ co/ \ CO/ 

where superscripts ', ", '", represent the first, second, third, etc., derivative with respect to z. 
Letting 

. riTTZi 

sin 



Gi ~ /nnV (AW « + ,W J^y 



ft) / \co / 



nrrzi 
cos 



^(, (z0 _ gi+l(2i) _Mta^l + . 



7177 L /7WT 

CO 

>i7r\ 2 /yi7r\ 4 
co / \co / 



j^ = gmM g&(q> 



(f r &)' ' 
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the 



2 2 "' ! 

A n =— {gi(Q)-B+(-l) n [Cco + B-g ni (w)\-H + (-l)»K}+ - V ft. 



For the condition that Z? = gi(0) and CcL)-hB=g m (co), the Fourier coefficients are 



A„ = —{(-l)"K-H}+-"y G,. 

77/7T ft) — — 



(3) 



Assume for the region, ^ z ^ co, that v as a function of z behaves as shown in figure 1, where 
gu gi, gz, and g4 are known functions off, (perhaps as determined from experiment) and/1,/2, and 
f.i encompass regions of uncertainty or regions which cannot be defined by experiment. If it is 
known that there are no discontinuities within the regions of the uncertainties, it is possible to 
smooth or bridge between the known functions with a cubic polynomial which gives continuity to 
two adjacent g functions and their first derivatives. Examining the first region that includes an 
uncertainty, the cubic polynomial, f\ (z) = M -f M\ z~\- M 2 z 2 -\- M : :i z s is solved for the coefficients, M, 
such that giUi)=/i(zi), gi(zi)=fi(zi), g 2 (z 2 )=fi{z 2 ) and #2(^2) =// (£2), from which the necessary co- 
efficients are 



Mo 



M-,=- 



3[gAz>) - gi(zi)]z - e[(3z + e)g[(z x ) + (Sz- €)g!(z 2 )\ 
4e 3 

e[g[(zi) + g! z (z 2 )] +gAzi)-g2(z 2 ) 

4e : * 



where z = Z\ + € = z 2 — 6 is the midpoint between Z\ and z 2 . Assuming bridging between g\ and g 2 , 
then 



G,=F 1 +F 2 



(4) 



where 



Fi = 



1 



nrr\ 4 
a)) 



{gl"( 



£(*) 



mry 
a) 



Jl7r\ 3 
CO 



gi'(Zl) 



. . . [ sin 

& (*■) 



+ . . . f COS 



(mry J to 



TITTZ\ 
CO 



«i'(ft)- 



£M 



Jin 

CO 



...} 



mrz 2 



CO 

(4a; 



Figure 1. Variation of v with respect to z showing the 
defined functions gi, £:>. g3, and gi connected by functions 
fi. f 2 , rtrcd fa through regions of uncertainty. 
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an 



d 



F, l 



nrr 

CO 



fi (zi) sin /i (z 2 ) sin 



CO CO 



+ l 



Jl7r\ 3 
CO 



rn/ N riTTZ2 , , v mrz\ 

J\ (z-?) cos Jl(Zi) COS 



CO CO 



Making the substitutions Zi = z — €, z 2 =2 + e, f" (z) = 2Af 2 + 6M. } z, and // " (z) = 6M 3 in F 2 , yields 

TT /mre\ . nifz 
U sin 

^ (n7re\ rnfz 
1 I Icos 

-ki(zi)-*(*) + efe{ (ft) +«*(*)}] -^^ — (4b) 

7177 



CO 

where 

"w-^-i-fj+f;-... 

7 T (iQ)=^(sin i 8-)8cosi8) 

= q M fi 2 | ff 4 

\3-l! 5-3! 7-5! * ' 

Plots of U(f3) and T(/3) are given in figure 2. 

For the interval zi ^ z ^ z 2 , (4b) is the bridging or smoothing function that gives continuity 
to g\{z) and g 2 (z) and their first derivatives at Z\ and z 2 , respectively. This function also eliminates 
computation of the cubic polynomial /i(z), where the values of M 2 and Ms sometimes give inordi- 
nately large values with subsequent loss of significance. The same analysis can be applied to the 
regions encompassed by/ 2 (z) and/s(z) so that G 2 and G 3 can be found. 

3. Specific Cases 

Where successive differentiation terminates, the functions g\(z) may be represented by poly- 
nomials. There are many variations of gi(z) possible, but probably the best example is to examine 
a set of straight lines as given in the following specific cases. 

3.1. Functions as Straight Lines 

The functions g\{z) are represented by a set of straight lines. Here g\ , {z) = H = K = 0, so that 
(3) becomes 

A n =-YG f . (5) 

CO ^ 
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FIGURE 2. I ariation of\ (/3) and T(/3) with p. 




Figure 2. Variation of\J(/3) and T(/3) with /3. 

If a set of straight lines which do not intersect at critical positions z/, as shown in figure 3 are 
assumed, then 



~ I 0) \ r ,. . # / nt . W7TZ; ft) r . N ., 717TZ/ 

k/ = — k/(2/)— g\+\(zi)\ sin [g-/(z/)~ g, +1 (z,)Jcos , 

\mrj 0) titt co 



(6) 



where the functions and derivatives are evaluated at z\. For the cases where the lines do intersect 
at z/, i.e., g;=gi + ], the second term on the right side of (6) is zero. 



3.2. Straight Line Functions Smoothed Near Intersections 

The convergence of (6) when substituted in (5) and (2) is not very satisfactory, and it is possible 
to increase the convergence by joining the straight lines with smoothing curves that give continuity 
to the lines and their first derivatives at points z\ — y\ and Zi + €u as shown in figure 4. 

Proceeding as in the derivation of (4), where F\ = 0, yields 



d = [(€/ - yd (gi ^(Zi)-g'i(zi)) + 2(g i+l (zi) -gi(zi))] 



T{nXj) cos a,- 



-{gi +1 (Zi)-g'l(Zi}} 



U{nkj) sin a, 
/rm\ 2 



(7) 



where 



n7T to _u \ 






As e/ = y, — » 0, (7) becomes (6). There are many possible specific cases for (7), i.e., 6j = 0, 
y* = 0, e,=7i, etc. Essentially, smoothing increases convergence by introducing the factors T(J3) 
and f/(j3), j8 = n\/, to the terms on the right side of (6). The plots T{fi) and C/(j8) (fig. 2) indicate that 
the advantage of the increase in convergence is not realized until /3 > 0. 
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FIGURE 3. Variation of v as a set of straight lines with 
function discontinuities at z\ and z-i. 



FIGURE 4. Variation of v as a set of straight lines with 
smoothing in the regions of L\ and z>. 



3.3. Example of Smoothing 

Figure 5 is used as an example of the variation of the dependent variable v with the dimension 
z in order to show the need for smoothing in the region of a discontinuity. For a hypothetical 
process, consider the potential on the surface of a system of two materials which is surrounded 
by a third material (not a perfect insulator), and where, on the surface of the system, the region 
^ z ^ 0.5 represents the potential in a perfect conductor, the region 0.5 ^ z ^ 1.0 represents that 
in a material of finite conductivity, and the discontinuity represents the contact resistance between 
the two materials. Usually, experimental processes allow measurement of potential in the two 
regions, but measurements at or near the discontinuity are not always possible or desirable. As- 
suming homogeneous isotropic materials, measurements taken in the two separate regions give 
linear relations which, when extrapolated to the point of contact (z = 0.5), do not meet in a common 
point and give an apparent discontinuity. When it is known that the process must be continuous, 
i.e., the surrounding material is continuous at z = 0.5, then a smoothing function may be applied. 

Numerical values for substitution in (7) are zi = 0.5, gi = 1.0, #1 = 0, #2 = 0.95, #2 = 1.9, and 
e = y, where a = /i7r/2, \ = tt€. Substitution in (7), (5) and (2) gives 



i-z-^(-ir 



T(2nk) sin lirnz , 3.8t/[(2n- 1)A] sin (2n- l)irz 



20^77 



(2rc-l) 2 7r 2 



(8) 



Figure 6 is a plot of v in the region of the discontinuity shown in figure 5 for the summation of 200 
terms of the above series and values of e = 0.0, 0.002, 0.003, and 0.004, and their respective smooth- 
ing curves. Figure 6 shows the influence of the increased convergence of the series toward its 




Figure 5. Variation of v with z for the problem of 
section 3.3. 
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Figure 6. Numerical values of v from (8) in the region 
of the discontinuity of figure 5 for a summation of 200 
terms and values ofe = 0, 0.002, 0.003 and 0.004. 
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intended values witli increase in the value of e. For this case, fairly good convergence of the series 
is attained for n > l/2e. For the summation of 300 terms, the fluctuations shown on figure 6 are 
virtually eliminated, except for e = 0. 

4. Discussion 

The example of section 3.3 as shown in figure 5 was chosen because it illustrates an extreme 
condition that a physical process may undergo; that is, the dependent variable has a zero slope 
up to the discontinuity at the point z = 0.5. It is then to be expected that actual physical processes 
where apparent discontinuities occur would give better convergence when smoothing conditions 
are employed. 

Generally, the Fourier sine series (2) can be used to represent a piecewise continuous function 
over an interval, ^ z ^ w, by the use of (3). Numerical evaluation will, of course, present prob- 
lems of convergence in the regions where discontinuities occur (as shown in fig. 6, e = 0). The 
lack of sufficient convergence in (2) can cause serious problems when (2) is used in the solution 
of (1). It is therefore expedient to employ the bridging or smoothing relationship (4) (or (7) for 
the special case of a set of straight lines) for which adequate convergence can be assured for a 
reasonable number of terms in the series of (2). 

A fairly simple heat transfer system that may be analyzed employing the method presented 
in this paper is shown in figure 7a by a cross section through a composite cylinder. The inner 
cylinder (regions 1 and 2) ^ r ^ a is composed of two materials of different thermal conductivities, 
usually with values greater than that for the material of the outer cylinder (region 3). For the 
determination of thermal conductivity, region 1 may comprise a material of known thermal con- 
ductivity and region 2 a material of thermal conductivity to be determined. It is necessary to 
determine the heat transfer characteristics in the material of region 3 in order to evaluate the radial 
heat loss or gain of the inner cylinder, ^ r ^ a. 

The temperature distribution at r = a will vary according to the thermal conductivities of the 
materials of regions 1 and 2, and the distribution may have an infinite number of variations. One 
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FIGURE 7. #ea£ transfer system. 

(a) Cross-section through composite cylinder. 

(b) A possible variation in temperature potential, v at r=o and r=b. 




such distribution for r=a with a thermal contact resistance at z = 0.5 is shown in figure 7b, where 
the intended heat flow is in the positive z-direction with a unit temperature potential at z = and a 
zero value at z=\. On r=b the temperature is assumed to decrease linearly in the above limits, 
considering region 3 only, a solution of (1) (assuming no (/> dependence) is 



v =l-z+Z 



^A n Fo(n7rr, mrb) sin mrz 



F ( )(rnra, mrb) 



(9) 



where F () (x, y) = Io(x)K {) (y) — Io(y)Ko{x) and /<> and K {) are modified Bessel functions of zero order and 
first and second kind, respectively. The substitution of r = a in (9) gives a form similar to (2), for 
which a solution for A n can be found from (7), the temperature distribution shown in figure 7b and 
a reasonable value of e that insures convergence of the resulting series. 

Two of the objectives of the NBS Heat Transfer Laboratory are to develop and analyze methods 
for the determination of thermal conductivity of materials which vary widely in thermal properties 
and to make measurements of thermal conductivity for government agencies. For both analysis 
and experiment, the smoothing or bridging relationships described are applicable. 

(Paper 71C2-248) 
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